arXiv: 1502.05801 v2 [astro-ph.SR] 2 Mar 2015 


Astronomy & Astrophysics manuscript no. lowmassfinal © ESO 2015 

March 3, 2015 


Improved angular momentum evolution model for solar-like stars 

II. Exploring the mass dependence 

F. Gallet 1 ’ 2 and J. Bouvier 1,2 


1 Univ. Grenoble Alpes, IPAG, F-38000 Grenoble, France 

2 CNRS, IPAG, F-38000 Grenoble, France 
e-mail: florian.galletSobs .ujf-grenoble . f r 

Received accepted - 


ABSTRACT 

Context. Understanding the physical processes that dictate the angular momentum evolution of solar-type stars from 
birth to maturity remains a challenge for stellar physics. 

Aims. We aim to account for the observed rotational evolution of low-mass stars over the age range from 1 Myr to 10 
Gyr. 

Methods. We developed angular momentum evolution models for 0.5 and 0.8 Mq stars. The parametric models in¬ 
clude a new wind braking law based on recent numerical simulations of magnetised stellar winds, specific dynamo and 
mass-loss rate prescriptions, as well as core-envelope decoupling. We compare model predictions to the distributions of 
rotational periods measured for low-mass stars belonging to star-forming regions and young open clusters. Furthermore, 
we explore the mass dependence of model parameters by comparing these new models to the solar-mass models we 
developed earlier. 

Results. Rotational evolution models are computed for slow, median, and fast rotators at each stellar mass. The models 
reproduce reasonably well the rotational behaviour of low-mass stars between 1 Myr and 8-10 Gyr, including pre-main 
sequence to zero-age main sequence spin up, prompt zero-age main sequence spin down, and early-main sequence conver¬ 
gence of the surface rotation rates. Fast rotators are found to have systematically shorter disk lifetimes than moderate 
and slow rotators, thus enabling dramatic pre-main sequence spin up. They also have shorter core-envelope coupling 
timescales, i.e., more uniform internal rotation. As for the mass dependence, lower mass stars require significantly 
longer core-envelope coupling timescales than solar-type stars, which results in strong differential rotation developing 
in the stellar interior on the early main sequence. Lower mass stars also require a weaker braking torque to account for 
their longer spin-down timescale on the early main sequence, while they ultimately converge towards lower rotational 
velocities than solar-type stars in the longer term because of their reduced moment of inertia. We also find evidence 
that the mass dependence of the wind braking efficiency may be related to a change in the magnetic topology in lower 
mass stars. 

Conclusions. We have included in parametric models the main physical processes that dictate the angular momentum 
evolution of low-mass stars. The models suggest that these processes are quite sensitive to both mass and instantaneous 
rotation rate. We have worked out and reported here the main trends of these mass and rotation dependencies, whose 
origin still have to be addressed through a detailed modelling of magnetised stellar winds, internal angular momentum 
transport processes, and protoplanetary disk dissipation mechanisms. 

Key words. Stars: evolution - Stars: solar-type - Stars: low-mass - Stars: rotation - Stars: mass-loss - Stars: magnetic 
field 


1. Introduction 

Observational constraints on the rotational evolution of 
low-mass stars have exponentially increased in the last 
decade, thanks to a number of dedicated long-term mon¬ 
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evolution. To account for the observations, parametric mod¬ 
els have to incorporate at least three major processes: the 
star-dis k interaction during the early pre-main sequence 
fPMS) ( Matt fc Pudritd 2005bl; IZanni fo Ferreiral 20091: 


Ferreira et al.1 120001 : iMatt et al.1 i201d : Zanni fc Ferreiral 


et al I2014f ). These new observational results provide ex¬ 
tremely useful guidance for the modelling of angular mo¬ 
mentum evolution of lo w-mass stars (M* < 1.2M© from 
1 Myr to 10 Gyr fe.g. fl rwin et all 120071 : iBouvierl 120081 : 

lDenissenkov_etjiI| j20jll|^JS2 adae 1^aljj201T|_Reiners & 

Mohantv l2012l : iGallet fc Bouvied 1201311 and offer a unique 
insight into the physical processes that dictate rotational 
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l201ll . 120121 : ICharbonnel et ahl 1 201 .'ll ) . The angular veloc¬ 
ity evolution of low-mass stars appears to be controlled to 
a large extent by these three physical processes, from their 
birth to the end of their MS evolution a nd possi bly beyond. 

In the first paper of this series dGallet fc Bouviej[2013fl , 
we developed a model to account for the angular momen¬ 
tum evolution of solar-nrass stars. The parametric model 
included star-disk locking, a new braking law with dynamo, 
and mass-loss pr escriptions that relied on the latest numer - 
ical simulations (ICranmer fc Saai1l20 1 ll : iMatt et, ahll20 12 ali , 
and core-envelope decoupling to account for the reduced ef¬ 
ficiency of internal angular momentum transport processes. 
The aim of the present study is to extend this model to 
lower mass stars: Can similar parametric models account 
for the evolution of lower mass stars? How do model pa¬ 
rameters depend on mass? What does this mass dependence 
reveal about the underlying physical processes? These are 
the questions we attempt to address in this study. In Sect. [2j 
we present the set of 18 rotational period distributions that 
we used to constrain our models at various ages. We briefly 
describe in Sect. [3] the assumptions we made to compute 
the angular momentum evolution of low-mass stars from 
birth to the end of the MS. The results are presented in 
Sect. [2 In Sect. [5j we work out the mass dependence of an¬ 
gular momentum evolution and discuss its implication for 
the physical mechanisms involved. We will also discuss the 
impact of these new models on gyrochronology and lithium 
depletion. Conclusions are drawn in Sect. [6] 

2. Observational constraints on the rotational 

evolution of low-mass stars: the data sets 

The goal of the angular momentum evolution models is 
to reproduce the observed rotational evolution of low-mass 
stars from birth to the end of the main sequence. To char¬ 
acterise the latter, we used rotational period distributions 
measured for the coeval populations of 18 star-forming re¬ 
gions and open clusters, thus spanning an age range from 
1 Myr to 1 Gyr, to wh ich we added results pe rtaining to 
old disk field stars from iMcQuillan et al.1 (12014 ). The data 
for star-forming regions and young clusters, including their 
age and angular velocity percentiles (25 th , 50 th , and 90 th ) 
in each mass bin investigated here (0.4-0.6, 0.7-0.9, and 0.9- 
1.1 M 0 ), are provi ded in T able HU see Appendix A of Gallet 
& Bouvier T2013l . for a detailed description of the properties 
of these clusters). 

Figure [T| shows the distributions of rotation periods 
plotted as a function of stellar mass for star-forming regions 
and young open clusters. Several patterns can be recognised 
in this figure that characterise the rotational evolution of 
young low-mass stars. The youngest star-forming regions, 
with an age less than 10-15 Myr, exhibit a wide distribution 
of rotation periods, ranging from less than 1 and about 10 
days, with little dependence on mass over the mass range 
investigated here (0.4-1.1 M 0 ). A well-defined relationship 
between rotation and mass starts to be seen in older clus¬ 
ters, from the ZAMS onwards, which first appears for solar- 
mass stars and later propagates to lower masses. The oldest 
clusters indeed exhibit quite a tight rotation-mass correla¬ 
tion, with rotation rates steadily decreasing towards lower 
masses. It is also noteworthy that the dispersion of rota¬ 
tional velocities is the lowest at ZAMS for very low-mass 
stars (cf. M50, 130 Myr) and increases later-on on the early 
MS over a few 100 Myr, while it is quite the opposite for 


solar-type stars whose wide rotational dispersion at ZAMS 
(cf. Pleiades, 125 Myr) is promptly erased on the early MS. 
This complex behaviour reflects the combination of various 
physical processes acting on the angular momentum con¬ 
tent of the stars as they evolve. 

Figure [2] shows the same data set as in Fig. |T] except 
that for each cluster, we extracted the observed angular 
velocity distribution in a given mass bin and plotted it as 
a function of time. The three panels of Fig. |T| thus illus¬ 
trate the evolution of the rotational distributions of 0.5, 
0.8, and 1 M 0 stars, respectively. The 25 th , 50 th , and 90 th 
percentiles of the distributions and their respective error 
bars were computed at each mass and age as described in 
iGallet fc Bouvierl (120131) . Figure [2] suggests a qualitatively 
similar evolution of the rotation rates for stars in the three 
mass bins. The evolution is characterised by a spin rate 
that hardly or slowly varies during the first few Myr, then 
followed by a rapid acceleration towards the ZAMS, and 
a subsequent braking over longer timescales on the MS. 
Quantitatively, the rotational evolution of the lowest mass 
stars appears delayed compared to solar-type stars, as they 
reach the ZAMS later and are braked on the main sequence 
over longer timescales. The models presented in the next 
sections aim to reproduce the overall evolution as well as 
subtle differences as a function of mass. 

While these distributions appear statistically robust and 
relatively unbiased, some residual contamination by non¬ 
cluster members cannot be totally dismissed. Indeed, some 
stars appear to rotate beyond break-up, presumably due 
to contamination of the samples by short period field bina- 
riet3. The reported photometric period is sometimes associ¬ 
ated with harmonic of the rotation period e.g. in case of two 
spots located at opposite longitudes on the stellar surfacefl 
Finally, some stars simply have a relatively low member¬ 
ship probability or an age estimate in conflict with that of 
the clustefl 

Another issue lies in the difficulty of estimating the 
age of clusters. Thus, for very young clusters with an age 
less than about 5-10 Myr, age uncert ainties may be ap - 
proximately the actual age estimate dBell et al.ll 120118 ). 
Fortunately, the change in rotation rates at these very 
young ages are relatively mild, so that the modelling does 
not suffer much from these uncertainties (i.e. shuffling the 
youngest clusters within the 1-5 Myr age range in Fig. [2] 
would not significantly impact the models). The ages of 
older clusters is much better known, with uncertainties 


1 We thus rejected M50-5-1624, -3-1468, -3-464, -7-7623, - 

3- 5840, -7-5624, -5-2673, -3-7334 -3-2531, -6-1574 -4-2077, - 

4- 4939, and -8-6076 from M50 (llrwin et all 120091) . N2362-2- 
6989, and N 2362-5-4947 from NGC 2362, #11041 (number from 
iHillenbrandll 19971) fr om ONC. and #494 (numbe r from Moraux 

et al. l2013l) from h Per, and #3245 (number from lLittlefair et all 

l2010i ) from CepOB3b 

2 We thus rejected 1SWASP J083722.23+201037.0 (KW 30) 
and 1SWAS P J0840 05.72+190130.7 (KW 256) in Praesepe 
( Delorme et al „ 201 1 ). 

3 We thus rejected JS 634 in Praesepe that has a mernber- 
ship prob a bility of 0 .62, compared to 0.94 for cluster members 
(lAgiieros et ^ 12011 . 

4 In NGC6530 (2 Myr), we rejected XID 138 (0.94 Mg, P ro t = 
0.18743 d) as its age is estimated to be 2 2 Myr a nd i s probab ly 
already accelerating towards the ZAMS dHenderson fc Stassunl 
120121 ) . 


2 














































Table 1 . Star-forming regions and open clusters whose rotational distributions are used is this study. 


Cluster 

Age 

(Myr) 

Ref. 

1 Mq 

H 25 

0.8 Mq 

0.5 Mq 

1 Mq 

^■50 

0.8 Mq 

0.5 Mq 

1 Mq 

H 90 

0.8 Mq 

0.5 Mq 

ONC 

1.5 

1 

3.28±0.3 

2.86±0.14 

3.59±0.26 

4.58±0.86 

4.03i0.58 

5.27±0.89 

13.21i2.14 

20.llil.93 

14.93i3.63 

NGC 6530 

2 

2 

3.87±0.55 

2.89±0.19 

3.19±0.21 

6.15i0.9 

4.08±1.14 

5.74±0.54 

22.32i5.15 

18.73i2.73 

22.72i3.48 

NGC 2264 

3 

3 

- 

2.81i0.5 

2.81i0.38 

- 

5.32i0.91 

6.16±1.93 

- 

10.66i2.01 

27.57i2.83 

Cep OB3b 

4 

4 

3.25±0.39 

3.91±0.81 

3.38±0.2 

6.99±1.07 

6.45±1 

5.58±0.43 

14.98il.55 

17.13i2.14 

23.66i5.42 

NGC 2362 

5 

5 

3.13±0.29 

2.75±0.29 

3.9±0.26 

4.2±0.58 

3.82±0.29 

6.31±0.67 

12.36i4.76 

14.01i4.4 

23.09i5.94 

h Per 

13 

6 

4.75±0.19 

3.93±0.17 

4.5±0.53 

8.31±0.85 

6.33±0.74 

11.93±5.03 

73.37i6.32 

46.99i6.2 

58.4i4.37 

NGC 2547 

35 

7 

- 

4.55±0.38 

4.77±0.87 

- 

5.28±1.01 

10.82±1.86 

- 

47.31i22.35 

45.24il7.7 

IC 2391 

50 

8 

- 

4.38i0.76 

- 

- 

9.76±4.29 

- 

- 

98.13i31.77 

- 

a Per 

80 

8 

- 

5.88±0.75 

- 

- 

14.81i8.32 

- 

- 

98.41il4.26 

- 

Pleiades 

125 

9 

4.92±0.13 

3.59±0.08 

7.98±1.2 

6.29±0.31 

4.27i0.24 

18.24±1.91 

37.99il0.08 

61.49i5.5 

76.74i4.91 

M 50 

130 

10 

3.38±0.3 

3.21±0.1 

6.79±0.53 

5.17±0.36 

4.23±0.22 

17.19±1.74 

15.41i5.16 

43.12i3.55 

72.22i3.32 

M 35 

150 

11 

4.39±0.07 

3.38±0.1 

- 

5.05±0.22 

4.59±0.71 

- 

24.3i5.27 

59.3i5.55 

- 

NGC 2516 

150 

12 

- 

- 

5.11±0.66 

- 

- 

11.52±1.56 

- 

- 

70.12i4.39 

M 34 

220 

13 

3.68±0.14 

2.67±0.1 

- 

4.64±0.83 

3.07i0.09 

- 

27.02i6.78 

21.31i5.14 

- 

M 37 

550 

14 

2.95±0.05 

2.38±0.04 

1.97±0.44 

3.26±0.06 

2.63±0.06 

Il.23i3.37 

3.75i0.13 

3.7i0.27 

40.69i2.14 

Praesepe 

580 

15 

2.73±0.05 

2.11±0.05 

1.39±0.07 

2.85±0.02 

2.52±0.11 

3.28il.57 

2.92i0.75 

5.94il.39 

22.03i4.63 

Hyades 

625 

16 

2.44±0.03 

2.05±0.04 

1.41±0.07 

2.62±0.06 

2.18i0.03 

1.91i0.09 

2.91i0.1 

2.46i0.15 

2.3i0.78 

NGC 6811 

1000 

17 

2.27±0.04 

- 

- 

2.36±0.01 

- 

- 

2.45i0.05 

- 

- 


et al. (l2013f) : (Ylllrwin et al.l (l2008bb : f81 Irwin et Bouvier; ('9)lHartman et all (120101): flOillrwin et al.l (l2009f) : ('ll ) iMeibom et all (120091); 1121 llrwin et al.l (l2007f) : f 131 " Meibom 

et al. <1201 lbfl : 1141 lHartman et all (1200911 : ('15) lAgiieros et al.l (1201 ill 4 lDelorme et all <1201 lh : 116) iDelorme et all 1120111 ): flT) IMeibom et, all (l2011ah .~ 
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Fig. 1. The rotational period distribution of 0.4-0.6, 0.7-0.9, and 0.9-1.1 M 0 stars (respectively shown as tilted triangles, 
crosses, and tilted squares) in star-forming regions and young open clusters is plotted as a function of stellar mass. The 
panels are ordered by increasing age, from top to bottom and left to right. The few black tilted crosses in NGC 6530, 
M50, and Praesepe are stars rejected from the samples for reasons explained in the text. The rotational periods of stars 
outside the selected mass bins are shown as small dots, as are subsamples not rich enough (i.e. N sta r < 15) to yield 
statistically meaningful percentiles and therefore rejected from the analysis. 


close to 10%. We therefore neglect the age uncertainties 
in the models developed below. 

3. Parametric models of angular momentum 
evolution: the assumptions 

The angular momentum evolution of isolated low-mass 
stars is controlled by the balance between three main phys¬ 
ical mechanisms: angular momentum removal by magne¬ 
tised stellar winds (here after “the wind braking”), the 
star-disk interaction, and the angular momentum transfer 
within the stellar interior. We discuss in this section the 
corresponding model assumption. 

3.1. Structural stellar evolution 

We adopt the iBaraffe et al.l (11998H non-rotating, solar- 
metallicity models computed for 0.5, 0.8, and 1 M 0 stars, 
with a mixing length parameter a=1.9 and an helium abun¬ 
dance Y = 0.275. These models provide the structural evo¬ 
lution (i.e. mass, radius, moment of inertia) of the inner 
radiative core and outer convective envelope from the early 
PMS to the end of the MS i.e. from 1 Myr to 10 Gyr. 

Figure [3] shows the evolution of the internal structure 
for 0.5, 0.8, and 1 M 0 stars. Most stellar quantities increase 


for higher mass stars (i.e. the stellar radius and moment of 
inertia, the moment of inertia of the radiative core, as well 
as its radius and mass), except for the moment of inertia of 
the convective envelope that decreases for higher masses. 
The size and mass of the radiative core is strongly mass 
dependent. While for a 1 M 0 star the radiative core rep¬ 
resents about 73% of the stellar size, this fraction falls to 
60% for a 0.5 M 0 star. Similarly, for a 1 M 0 star almost all 
the stellar mass is contained in the radiative core while only 
73% of the total mass is stored in the core of a 0.5 M 0 star. 
Note that the rapid decrease of the star’s moment of iner¬ 
tia during PMS evolution is expected to strongly impact 
on the star’s surface velocity, since angular momentum is 
given by J = I ■ 12 where I is the moment of inertia and f 1 
the angular velocity. 


3.2. Magnetised stellar winds 

Assuming a spherical outflow, the braking torque exerted 
by a magnetised stellar wind on the stellar surface can be 
expressed as 


^ wind OC * M w i nc [ • (1) 
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Fig. 2. Angular velocity distributions are plotted as a func¬ 
tion of time for low-mass stars in star-forming regions and 
young open clusters. Each panel covers a different mass 
bin: 0.9-1.1 M 0 (upper panel), 0.7-0.9 M 0 (middle panel), 
0.4-0.6 M 0 (lower panel). The red, green, and blue tilted 
squares and associated error bars represent the 25 th , 50 th , 
and 90 th rotational percentiles, respectively. The open cir¬ 
cle shows the angular velocity of the present Sun for refer¬ 
ence. The black rectangle labelled OD (lower right corner 
of each panel) shows the angular velocity dispersion of old 
disk field stars. The black dashed line represents the evolu¬ 
tion of break-up velocity. The left vertical axis is labelled 
with angular velocities normalised to the Sun’s, while the 
right vertical axis is labelled with rotational periods (days). 


where M W i nd is the mass-loss rate, and xa is the averaged 
value i rftheA]^eiy^ l dms fii The - latter_is_obtainedJrOTn i Matt 
et al. (1201281 1 as 


r A = Ki 


n 2 td2 

p**'* 


M wind y/K'ivl c + niR'i 


lit 


( 2 ) 








Fig. 3. Evo lution of the main stellar parameters from the 
iBaraffe et al.l (ll998t )’s models for 1 M 0 (solid line), 0.8 M 0 
(dotted line), and 0.5 M 0 (dashed line). Upper panel, from 
left to right: stellar radius (in solar unit), moment of inertia 
(in solar unit, 7 0 = 6.411 x 10 53 g.cTO -2 ), and moment of 
inertia of the convective envelope (in stellar unit). Lower 
panel, from left to right: the radiative core’s moment of 
inertia, radius, and mass (in stellar unit). 


where B p is the surface strength of the dipole magnetic field 
at the stellar equator, v esc = \/2GM* /R * is the escape ve¬ 
locity, and m = 0.22, K 2 — 0.0506. As in lGallet fc Bouvieil 
(120131 ). we identify B p to the strength of the mean magnetic 
field B„f* where B * is the magnetic field intensity and /* is 
the filling factor, i.e. the fraction of the stellar surface that 
is magnetised. The evolution of the mean magnetic field 
B*f * as a function of stellar density, effective temperature, 
and an gular vel o city is taken from the model developed by 
ICranmer fc Saaii (1201 lh . To reproduce the mean solar filling 
factor (/ 0 = 10” 1 — 10~ 3 ), we s lightly modified the expres¬ 
sion of ft given by ICranmer fc Sa.a.il (1201111 , as described in 


iGallet fc Bouvier (120131 ). Eventually, this yields 


B*f* = 1.13 x B eq x /*, 

/r 

= 1.13a 


l8nptk B T e ff 


0.55 


A im H [1 + (x/0.16) 2 - 3 ] 


1.22 ’ 


( 3 ) 


where x = Ro/Roq and Ro = P ro t/T CO nv is the Rossby 
number, i.e. the ratio of P ro t the rotation period to r conv the 
conv ective turnover timescale , with Roo = 1.96 ( Wright 
et al. 120111 : iQglethorpe fc Garaudl 12013^ . In Eq. [3l B eq 
is the_egui£artilion_jnagnetic_^eld_strength_{Cranmer & 
Saar 1201 il l with p* the photospheric mass density, k B the 
Boltzmann’s constant, T e ff the effective temperature, p the 
mean atomic weight, and to h the mass of a hydrogen atom. 

The ma ss-loss rate prescription we inject in Eq. [5] also 
stems from ICranmer fc Saaii (120111 ) coronal wind models 


M w ind OC 




-2/7 

(Ft.) 12 '" 


.5/7 

J* 5 


( 4 ) 


where A* is the stellar luminosity, and Fa* is the quantity of 
energy deposited in the photosphere by Alfven waves. This 
energy is converted into a heat energy flux that powers the 
stellar wind. We also modified the original expression of Fa* 
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from lCranmer fc Saail ( 2011 ') to incorporate the difference 
of mixing len gth parameters used in t heir study and ours 
(see details in Gallet fe Bouvieil (l2013f )l. 

With these assumptions, Fig. [4] shows the evolution of 
R*/*, M W i n d, and J as a function of the angular velocity. 
All three quantities exhibit different regimes, with a steeper 
dependence on velocity at low spin rates and a shallower 
slope at high velocities. Indeed, this reflects the saturation 
of /* at high velocities (see Eq. [3]). At a given velocity, 
lower mass stars exhibit higher mean magnetic fields, which 
primarily results from their larger photospheric pressure. 
The magnetic field also appears to saturate faster in lower 
mass stars. This is a consequence of low-mass stars hav¬ 
ing longer turnover timescale, while saturation occurs at 
the same Rossby number for all masses ( Ro sa t — 0.1, cf. 
Eq. [3] and the inset in Fig. [JJsjipper panel). Unlike previ¬ 
ous studies that followed [Kawaleil f l988 Ts prescription and 
used distinct relationships to account for the non-saturated 
and saturated magnetic regimes, we use a single prescrip¬ 
tion that smoothly bridges the two regimes. Moreover, as 
this prescription scales on the Rossby number, it can be 
straightforwardly extrapolated to any mass range. 


Like the magnetic field, and for the exact same reason, 
the mass-loss rate first increases with angular velocity at 
low spin rates and then saturates at higher velocities. The 
strong mass dependence of the mass-loss rate seen in Fig. [4] 
stems from Fa* being very sensitive to mass. Fa* is a func¬ 
tion of T e ff and log<? only and rapidly decreases towards 
lower masses. The angular momentum loss rate can thus 
be computed at any age step for any mass by combining 
Eqs. [U El [21 and[H together with the structural evolution 
models of lBaraffe et al.l (ll998l h Figure [4] shows the run 
of angular momentum loss rate as a function of angular 
velocity. The same trend for a steep increase with rota¬ 
tion at low spin rates followed by a shallower rotational 
scaling at higher velocities is seen, which naturally results 
from the combined rotation dependencies of the magnetic 
field and mass-loss rate discussed above. The asymptotic 
express ions of the braking law are provided in Gallet & 
Bouvier ( 20131 ). Quantitatively, however, we stress that we 
had to renormalise the I\\ constant appearing in Eq. [2] in a 
mass-dependent manner to reproduce the observations (see 
Section 15.31) . With this renormalization, shown in Fig. EJ 
the braking torque is found to be roughly the same for all 
masses in the non-saturated regime, but saturates faster 
and at lower values for lower mass stars. We emphasise 
that this behaviour of the braking torque indeed holds the 
keys of the rotational evolution of aging low-mass stars: it 
encodes both the longer spin down timescale of lower mass 
stars (due to the weaker braking torque in the saturated 
regime) and their lower final velocities (due to a higher 
ratio of braking torque to stellar moment of inertia in the 
saturated regime). We return to this crucial point in Section 


3.3. Star-disk interaction 

As apparent in Fig.[2j observations suggest that during the 
first few Myr of the PMS phase a physical process acts to 
prevent the stellar surface from spinning up. This behaviour 
is believed to result from the magnetic interaction between 
the young stellar object and its accretion disk, even though 
the actual physical process is not totally elucidated yet (see 


Matt fc Pudritzj 2005al. 2008 al ibi: IZanni fc Ferreiral 1201 ll . 
20131: Bouvierel^a^ 2013h 


As in lGallet &; Bouvieil (l2013f ) the model developed here 
assumes that the angular velocity of the stellar surface is 
held constant as long as the star accretes from its circum- 
stellar disk. The disk lifetime T<u s k , a free parameter of the 
model, thus dictates the duration over which the star is 
kept at its initial spin rate. When the disk eventually dis¬ 
sipates, after a few Myr, the disk coupling ends, and from 
thereon, angular momentum is removed at the stellar sur¬ 
face by stellar winds. 


3.4. Core-envelope decoupling 

The transport of angular momentum in the stellar interior 
is one of the most importa nt physical processes tha t occurs 
during the early MS phase. fGallet fc Bouvierl (I2013H showed 
that this process is central to the rotational evolution of 
solar-mass stars as it leads to the storage of a large fraction 
of the stellar angular moment within the rapidly spinning 
radiative core. As many uncertainties remain in the actual 
physical process tha t contro l ang ular mo mentum transport 
in stellar interiors dCharbonnel fc Talon 120051: Charbonnel 
et al. 120131b we adopt here a two-zone model where the in¬ 
ner radiative core and the outer convective envelope, each in 
uniform ro tation are loosely coupled. We follow MacGregor 
& Brenner (j 199 il l in assuming that a quantity A J is ex¬ 
changed between the core and the envelope over a timescale 
r c — e (hereafter the core-envelope coupling timescale). This 
quantity A J is the amount of angular momentum that the 
core and the envelope have to exchange to restore uniform 
rotation throughout the star, and is given by 


A J = 


Jp.nrp 3r.nrp.Jt 


L env° core 


core 0 env 


(5) 


where I and J refer to the moment of inertia and angular 
momentum of the radiative core and t he convective enve¬ 
lope, respectively. We also assume, as in lAllainl (|l998l h that 
T c — e is constant for a given model (see Sect. 3 of Gallet & 
Bouvier T2013i . for a more detailed discussion about r c _ e ). 


3.5. Equation of evolution 

The angular momentum evolution of low-mass stars is con¬ 
trolled by the internal and external torques applied on and 
within the stars (see Sect. El). This evolution can be ex¬ 
pressed as 


d J dl* dfl* 

dt dt * ”*" * dt 


( 6 ) 


where J is the angular momentum, / the moment of inertia, 
f l* the angular velocity, and T the sum of all the external 
torques. In Eq.[6]the quantity 7T2* is the contraction torque. 
In the framework of the assumptions described above, the 
equations that control the angular velocity evolution of the 
envelope can be expressed using Eq. [G] as 

if t A Tdisk- 


1 1 prim: — - I > 


if t > T di , 


sk • 


dfl r 


dJ ci 


dirt 


dt 


dt Ir. 


dt 


(7) 
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Fig. 4. Upper panel: Mean magnetic field strength as a function of stellar angular velocity for 1, 0.8, and 0.5 M 0 stars 
(solid, dotted, and dashed lines, respectively). The insert shows the mean magnetic field as a function of the Rossby 
number. The Sun’s range of £?*/* = 2—7.7 G is shown as a vertical bar. Middle panel: Mass-loss rate as a function of stellar 
angular velocity for 1, 0.8, and 0.5 M 0 stars (solid, dotted, and dashed lines, respectively). The Sun’s present day mass- 
loss is shown as an open circle. Lower panel: Angular momentum loss rate normalised to the Sun’s as a function of angular 
velocity for 1, 0.8, and 0.5 M 0 stars (solid, dotted , and dashed lines, respectively). We used J 0 = 7.169 x 10 30 g.cm 2 .s -2 , 
and S2 0 = 2.87 x 10 -6 s _1 dWeber fc David[l967 ll. 


where 


dJ f,i 


dt 


= -r 


wind 


+ r ce — r 


radi 


( 8 ) 


with T W i n d the wind braking torque (see Eq. [T]), T ce = 
A J/r c _ e the quantity of angular momentum transported 
from the core to the envelope over the timescale r c _ e , and 
T ra d the rate of angular momentum lost from the convective 
envelope during the formation of the radiative core. As the 
core develops, a fraction of the convective envelope becomes 
radiative and the rate of angular momentum transferred 
from the envelope to the expanding radiative core is 



conv 


dM ra d 

dt 


( 9 ) 


where dM ra d is the quantity of material that is con tained 
in a thin shell at a radius R ra d inside the star (lAllainlll998l) . 
Finally, the total torque applied on the stellar surface is 


dJ. 


A J 


j, 1 ‘wind T RradS^conv 

at T c —e O 


dM ra d 
dt 1 


( 10 ) 


Similarly, the angular velocity evolution of the core is com¬ 
puted as: 

dTlrad _ 1 A J | 2 Rrad 0 dM ra d 

li T ' o T ^ ^conv t * 

aZ 1-rad 7~c—e ^ J-rad aZ 

dirad ^rad ^ 

TTi T ’ 

ai -Lrad 

4. Results: confronting models to observations 

The free parameters of the model are adjusted so as to best 
reproduce the observations. These are: the initial rotation 
period at 1 Myr, Pi n it, the core-envelope coupling timescale, 
T c _ e , the disk lifetime, Tdisk > and the calibration constant 
of the wind braking law, K\. The best parameters found 
for the slow, median, and fast rotator models are listed in 
Table [21 

4.1. The rotational evolution of low-mass stars 

Figure [5] summarises the angular velocity evolution of slow, 
median, and fast rotators in the three selected mass bins, 
centred on 0.5, 0.8, and 1 M 0 . 


and the associated angular velocity evolution is 


4.1.1. Solar-mass stars 


dfl r 


dt 


1 A J 2 R: 


3 Ir, 


ad 


L conv 1 c—e 


a. 


dl c , 


dMrad 

’ dt 

r wind 


dt 


Solar- mass models were already presented in Gallet & 
Bouvier (l2013fl . The data set has been updated for this 
study by adding the results recently published for five ad¬ 
ditional open clusters, namely: Cep OB3b (4 Myr), IC2391 
( n ) (50 Myr), a Per (85 Myr), NGC2516 (150 Myr), and M34 
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Age (Myr) Age (Myr) 

(a) 1 Mq (b) 0.8 Mq 



Age (Myr) 

(c) 0.5 Mq 


Fig. 5. The angular velocity of the convective envelope (solid lines) and of the radiative core (dashed lines ) is shown as a 
function of time between 1 Myr and 10 Gyr for slow (red), median (green), and fast (blue) rotator models in three mass 
bins centred on 0.5 Mq (bottom panel), 0.8 Mq (upper right panel) and 1 Mq (upper left panel). The left vertical axis 
is labelled with angular velocity normalised to the Sun’s, while the right vertical axis is labelled with rotational periods 
(days). The red, green, and blue tilted squares and associated error bars represent the 25 th , 50 th , and 90 th percentiles of 
the observed rotational distributions at each sampled age. The black rectangle labelled OD (lower right corner of each 
panel) shows the angular velocity dispersion of old disk field stars. The open circle is the angular velocity of the present 
Sun shown for reference, and the dashed black line illustrates Skumanich’s (1972) relationship, 12 oc i -1 / 2 . 


(220 Myr), although not all of them have enough stars in a 
given mass bin to constrain the models (see Table [T]). We 
also restricted the mass bin to stars strictly within the 0.9- 
1.1 Mq mass range. As a result, the values of the model 
parameters have changed slightly, and are listed in Table El 
Compared to the values listed in I Gallet fe Don vied (1201311 , 
we fixed the Ki parameter to 1.7 for the three rotator mod¬ 
els (instead of 1.8 for the slow and median rotators, and 
1.7 for the fast rotators), and used a longer disk’s lifetime 
for the slow rotator model (9 Myr instead of 5 Myr). Other 
parameters did not change by more than 20% from the pre¬ 
vious study. These updates do not significantly affect the 
overall angular velocity evolution of solar-like stars and the 
results and discussion provided in iGallet fc Bouvieil (|2013H 
remain valid. 


4.1.2. Stars of 0.8 Mq 

Figure [5 (b) | and 5(a) show a qualitatively similar angular 
momentum evolution for 0.8 and 1.0 Mq stars. The initial 


rotation rates, 1.4, 6, and 9 days for rapid, median, and 
slow rotators, respectively, are fixed by the percentiles of 
the rotational distributions of the youngest clusters. The 
rapid spin up between the early PMS and the 13 Myr h 
Per cluster requires a disk lifetime as short as 3 Myr for 
rapid rotators, increasing to 7 Myr for slow rotators. To 
rapidly accelerate the convective envelope up to the ZAMS 
where Q,* AMS ~ 50 — 60 12q, rapid rotators must maintain 
nearly solid-body rotation, which requires a core-envelope 
coupling timescale as short as 15 Myr. In contrast, median 
and slow rotators have moderate velocities on the ZAMS, 
which calls for a longer core-envelope coupling timescale 
of 80 Myr, so that the envelope is readily braked while the 
core still spins fast. In these models, the rotational gradient 
between the core and the envelope reache s a m aximum of 
Af2/f2 ~ 2 — 3 at about 100 Myr (cf. Fig. |6(b)[ ). 

As in the case of a 1 Mq model, we find that the core¬ 
envelope coupling timescale is similar in slow and moderate 
rotators, and much longer than for fast rotators. In turn, 
such a long coupling timescale implies that the resurfacing 
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Age (Myr) 

(c) 0.5 M e 


Fig. 6. Upper panel: Velocity shear at the base of the convective zone (fl core — Ll env )/Ll env in the case of fast (blue), 
median (green), and slow (red) rotator models. Lower panel: Spin-down timescale (| J/J\) expressed in Myr, in the case 
of fast (blue), median (green), and slow (red) rotator models, a) 1 Mq, b) 0.8 M 0 , and c) 0.5 Mq. 


Table 2. Model parameters. 


M* — 1 Mq 


Parameter 

Slow 

Median 

Fast 

Pinit (days) 

8 

5 

1.4 

T c _e (Myr) 

30 

28 

10 

Tdisk (Myr) 

9 

6 

2 

Ah 

1.7 

1.7 

1.7 

M t 

, = 0.8 Mq 


Parameter 

Slow 

Median 

Fast 

Pinit (days) 

9 

6 

1.4 

r c _e (Myr) 

80 

80 

15 

Tdisk (Myr) 

7 

5 

3 

K\ 

3 

3 

3 

M, 

, = 0.5 Mq 


Parameter 

Slow 

Median 

Fast 

Pinit (days) 

8 

4.5 

1.2 

r c _ e (Myr) 

500 

300 

150 

Tdisk (Myr) 

6 

3.5 

2.5 

K\ 

8.5 

8.5 

8.5 


of the angular momentum stored in the radiative core will 
occur over an extended period. Indeed, this long-term out¬ 


ward transport of angular momentum partly compensates 
for the angular momentum losses at the stellar surface due 
to stellar winds. This results in a relatively shallow velocity 
decrease on the early main sequence, as observed for slow 
and moderate rotators that exhibit a clear change in the 

ng. 

easily 


spin-down slope between 100 and 500 Myr (cf. Fig 
In contrast, rapid rotators are spun down much more 
on the early MS, as the braking torque scales with surface 
velocity. 


The angular velocity of the slow, median, and fast ro¬ 
tators eventually converge towards the same spin rate at 
about 1 Gyr. By the age of the Sun, all the models pre¬ 
sented here exhibit a nearly solid body rotation and have a 
surface velocity lo wer th an the ; Sun’s, in agreement with the 
Kepler results (see lMcQuillan et al.ll2013l) . The lower termi¬ 
nal velocity of 0.8 M 0 models compared to the solar-mass 
models directly stems from the calibration of the braking 
law. As shown in Table 0 the K\ constant of the brak¬ 
ing law was increased from 1.7 for solar-mass stars to 3 for 
0.8 M 0 stars. Indeed, regardless of the detailed rotational 
history, the terminal velocity on the main sequence only de¬ 
pends upon the normalization of the braking law, thus de¬ 
creasing as the K\ constant increases. Note that with this 
renormalization, the braking torque of lower mass stars in 
the saturated regime remains weaker than for solar-mass 
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stars (cf. Fig. U), which simultaneously accounts for their 
longer spin-down timescale on the main sequence. 


4.1.3. Stars of 0.5 


The main trends of the rotational evolution of solar-mass 
stars can still be recognised in that of 0.5 M 0 stars, as 
shown in Fig.5. A clear difference, however, is that the evo¬ 
lutionary timescales become much longer at lower masses. 
In particular, the rotational spin down occurs over several 
Gyr on the main sequence for 0.5 M 0 stars, compared to 
only a few 100 Myr for solar-mass stars. 

The early PMS evolution is not very different from that 
of more massive stars: similar disk lifetimes, increasing from 
2.5 Myr for fast rotators to 6 Myr for slow rotators, are re¬ 
quired to prevent the lower mass stars from spinning up. 
Once the disks are dissipated, the spin up occurs on a longer 
timescale to the ZAMS, which is reached at about 120 Myr. 
As for higher mass stars, even though stellar winds are ac¬ 
tive during this phase, they are unable to balance the accel¬ 
eration torque due to decreasing stellar moment of inertia. 

In comparison to the 1 M 0 , the core-envelope coupling 
timescale increases by a factor of 10 or more in 0.5 M 0 
models, amounting to 150 Myr for fast rotators and 500 
Myr for slow and moderate rotators. One recovers the in¬ 
creasing coupling timescale from slow to fast rotators, as 
was the case for 0.8 and 1 M 0 models. This allows ini¬ 
tially fast rotators to reach high velocities at the ZAMS, 
while the convective envelope of slow and moderate rota¬ 
tors is already significantly braked before the star reaches 
the ZAMS (see Fig.5). Indeed, the increased core-envelope 
decoupling in slower rotators, coupled to their systemati¬ 
cally longer disk lifetime, is at the origin of the enhanced 
dispersion of rotational velocities on the ZAMS. 

In these models, the largest amount of differential rota¬ 
tion between the inner radiative core and the outer convec¬ 
tive envelope is reached at 500 Myr f or slow rotators and 
amounts to Afi/fi ~ 3 (cf. Fig. |6(c)[ ). The resurfacing of 
angular momentum being transported from the core to the 
enve lope on a timescale of <1 Gyr on the MS is still seen in 
Fig. 5(c) as a kink in the spin-down rate of slow rotators, 
even though it is not as pronou nced a s in higher mass stars. 
This is also highlighted in Fig. 6(c) ( lower panel) where a 
small plateau can be seen at « 400 Myr in the spin-down 
timescale of slow rotators. 

The angular velocity of the slow, median, and fast rota¬ 
tors models eventually converges towards the same rotation 
rate at a bout 3 G yr in a state of nearly uniform rotation 
(see Fig. |5(c)| and 6(c) I. As explained above for 0.8 M 0 
stars, the rotational convergence is delayed for lower mass 
stars due to the weaker braking torque in the saturated 
regime (cf. Fig. U). However, when rotational convergence 
is completed, lower mass stars exhibit lower terminal ve¬ 
locities, due to their reduced moment of inertia. Indeed, at 
low velocities the ratio of braking torque to moment of in¬ 
ertia increases for lower mass stars (cf. Fig. |4|), i.e. lower 
mass stars are more effectively braked in the non-saturated 
regime, which directly results in lower terminal velocities. 


4.2. Differential rotation 


The radial differential rotation rate is defined here by 
ADl/DL conv where AH = D core — fl env . This quantity mea¬ 


sures the velocity shear at the tachocline between the ra¬ 
diative core and the convective zone of the two-zone model. 
This is obviously a simplification of the actual rotation pro¬ 
file expected to arise in the stellar interior from angular 
momentum redistribution. Nevertheless, it provides a mea¬ 
sure of how much angular momentum may be stored in the 
stellar interior at any given time of its evolution. 

Figure [6] shows the evolution of differential rotation 
for the models described above. A few systematic trends 
emerge. Somewhat unexpectedly, slow and moderate rota¬ 
tors have systematically larger differential rotation rates 
than fast rotators. This is the direct result of the assump¬ 
tion of longer core-envelope coupling timescale for slower 
rotators (cf. Table [2]) . This assumption is required to ac¬ 
count for the widening of the rotational velocity dispersion 
from the early PMS to the ZAMS: initially fast rotators 
have to be spun up quite effectively while the surface of 
initially slow rotators must be spun down before they reach 
the ZAMS. The only way to achieve this, supposing realistic 
disk lifetimes (< 10 Myr), is to have a weaker core-envelope 
coupling in slow rotators than in fast rotators. The direct 
consequence of this assumption is the stronger differential 
rotation in slow and moderate rotators on the ZAMS, with 
AQ,/Dl conv ~ 2 — 3 compared to <1 for fast rotators. 

Another clear trend is that the peak of differential ro¬ 
tation rate occurs earlier in slow and moderate rotators 
than in fast rotators. While the former exhibit a maximum 
tachocline shear at or close to the ZAMS, the maximum 
occurs later on the main sequence for the latter (cf. Fig.6). 
This differential effect presumably results from the com¬ 
plex interplay between the spin-down rate, which is steeper 
for fast rotators, and the timescale for angular momentum 
distribution, which is shorter for fast rotators (shorter core¬ 
envelope coupling timescales, cf. Table [2]). Clearly, the evo¬ 
lution of differential rotation in stellar interiors is relatively 
complex and depends on the details of the models. 

Overall, the two-zone models suggest that a large quan¬ 
tity of angular momentum can be hidden in the stellar inte¬ 
rior over duration of a few 100 Myr and up to about 1 Gyr 
on the main seq uence. This has obvious consequences on 
stellar evolution (ICharbonnel et al]f20 13jh_the depletion of 
light elements such as Lithium t: Bouvie r ‘2 0081: Eggenberger 
et al. l2012al f. and the type of magnetic dynamos that may 
be ins trumental in young solar-type stars fVi dotto et al.l 
l2014af ). We return to these aspects in the next section. 


5. Discussion 

Thanks to the measurement of thousands of rotation pe¬ 
riods in star-forming regions and young open clusters, we 
are now able to accurately trace the angular momentum 
evolution of low-mass stars from the birthline (< 1 Myr) to 
the age of the Sun. Over the mass range considered here, 
three main phases of rotational evolution can be identi¬ 
fied: a nearly constant surface rotation rate during the first 
few million years of the PMS phase, a sharp increase to¬ 
wards the ZAMS, and a steady decline on the early MS 
on a timescale of a few 100 Myr. Furthermore, observa¬ 
tions suggest a large initial spread of rotation periods for 
the youngest star-forming regions. At the start of the PMS 
phase, rotational distributions exhibit periods ranging from 
1-3 days to 8-10 days. This initial spread increases further at 
ZAMS, with periods ranging from 0.2-0.4 days to about 6-8 
days. Eventually, the spread is erased on the early MS over 
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a timescale of a few 100 Myr, as surface rotation converges 
to a low, mass-dependent terminal velocity. The models de¬ 
veloped in the previous section show that these trends can 
be well reproduced with a small number of assumptions: i) 
a magnetic star-disk interaction that is believed to prevent 
the stars from spinning up during the early PMS, which 
results in a nearly constant rotation rate during the disk 
accretion phase; ii) angular momentum losses due to mag¬ 
netised stellar winds, a process whose magnitude depends 
on stellar mass, magnetic field, and rotation rate, and which 
dominates the rotational evolution of low-mass stars past 
the ZAMS; and iii) angular momentum redistribution in the 
stellar interior, which allows a large fraction of the initial 
angular momentum to be temporarily stored in the inner 
radiative core on the early MS. 


The combination of star-disk interaction, wind braking, 
and core-envelope decoupling fully dictates the surface evo¬ 
lution of low-mass stars. In the following sections, we dis¬ 
cuss the relevance of each model parameters and their de¬ 
pendence on mass. 


5.1. Initial conditions and disk lifetimes 


In the early PMS, the central model parameter is the disk 
lifetime, Tdisk- Since a constant rotation rate is assumed 
during the accretion phase, the longer the disk lifetime, the 
lower the rotation rate at ZAMS. Thus, disk lifetimes of a 
few Myr are required by the models to reproduce the wide 
rotational distributions observed at ZAMS. Recent mea¬ 
surements of infrared excesses and disk accretion rates in¬ 
dicate that most stars are born with a disk, while only half 
of them still have an accreting disk at a bout 3 Myr, and 
very f ew indeed by an age of 10 Myr (e.g. Hernandez et al.1 
l2f)08l : IWvattl 120081 : IWilliams fc Cieza 1201 if ). Observations 
thus suggest an average disk lifetime of 3-5 Myr , and a max¬ 
imum lifetime of order of 10-20 Myr (jBell et al.ll2013f ). This 
is quite consistent with the disk lifetimes required by the 
models, which range from 2 to 9 Myr (cf. Table [2]). 


The combination of the initial period, Pinit , and the 
disk lifetime, Tdi S k , primarily drives the surface velocity at 
ZAMS. There is some degeneracy in the models between 
these two parameters, as a longer initial period could be 
compensated by a shorter disk lifetime in order to yield the 
same velocity at ZAMS. To solve this degeneracy, we fixed 
the initial period based on the observed rotational distri¬ 
butions of the youngest PMS clusters. As a result, we find 
that disk lifetimes do not vary much with stellar mass but 
are systematically longer for slow rotators than for fast ro¬ 
tators in e ach mas s bin fcf. Table 121). As discussed in Gallet 
& Bouvier (120131 1 for solar-mass stars, this correlation be¬ 
tween disk lifetime and initial period, which is shown here 
to extend to lower mass stars as well, may be a sequel to 
the star-disk interaction process operating in the embed¬ 
ded phase. In other words, more massive disks would yield 
lower initial velocities and live longer. Hence, a distribution 
of protostellar disk masses might actually be the source of 
the dispersion of initial rotation rates of low-mass stars at 
the start of their PMS evolution. 


5.2. Core-envelope decoupling and the timescale for angular 
momentum transport 


Past the disk regulation phase, the shape of the gyrotracks 
is mostly dictated by the equilibrium between the wind 
braking torque and the internal redistribution of the angu¬ 
lar momentum. The dominant angular momentum trans¬ 
port process has not been securely identified yet, but sev¬ 
eral candidates have been pr oposed, namely hydrodynami- 
cal instabilities (see e.g. lKrishnamurthi et al.lll997ll. internal 
magnetic fields (see e.g. Denisscnkov & Pinsonneaultil2007l : 
IDenissenkov et all 12010 ; [Sj3ada^^_ah :2010j^Eggenberger 
et al. 2012af). and internal gravity waves (see Talon & 
Charbonnel l2005l 120031 : ICharbonnel et al.ll2013l f. 

Here, we used a simplified two-zone model consisting of 
a radiative core and a convective envelope, which are both 
in uniform rotation but at a different rate. The assumption 
of core-envelope decoupling thus leads to a velocity shear, 
indeed a discontinuity, at the tachocline. This obviously 
should be considered as a simplified approxi matio n of more 
complex internal rotational profiles (e. g. Snada et al.ll201f)l : 
IDenissenkov et all 120101: iBrun et, al.ll201ll: Turck-Chieze 
et al. 1201 It Lagarde et all 1 201 2t ) . As described in Sect. [3J 
we consider that the angular momentum transport mecha¬ 
nism acts to restore uniform rotation between the core and 
the envelope. The central parameter is the core-envelope 
coupling timescale, r c _ e , which is the time it takes to reach 
a state of uniform rotation. The efficiency of the angular 
momentum transport process is then measured by r~_ e . 

A striking result of our models is that r c _ e appears to 
strongly depend both on mass and on spin rate. As Table [2] 
summarises, lower mass stars have coupling timescales at 
least 10 times longer than solar-mass stars, and slow ro¬ 
tators have coupling timescales 3 to 6 times longer than 
fast rotators. Indeed, the varying coupling timescales are 
mostly responsible for producing the noticeable differences 
exhibited by the slow and fast rotator models presented in 
the previous sections. 


A small value of r c _ e for fast rotators will lead to 
high rotation rate at ZAMS since angular momentum is 
effectively transported from the core to the envelope over 
a timescale shorter than the contraction timescale. Once 
on ZAMS, rapid spin down will occur because of the f l- 
dependence of the wind braking law. This accounts for the 
triangular shape of the fast gyrotracks relative to the ZAMS 
(cf. Fig.5). In contrast, the much longer r c _ e of slow rota¬ 
tors results in a more complicated gyrotrack shape, which 
peaks slightly before the ZAMS as the envelope is spun 
down first and flattens out on the early MS as angular mo¬ 
mentum hidden in the core resurfaces in the envelope over 
several r c _ e . The main differences between slow/moderate 
and fast rotator models thus primarily arise from their 
widely different coupling timescales. It is worth noticing 
that this behaviour is true over the full mass range investi¬ 
gated here. 

Bouvier (l2008tl and llrwin fc Bouvier! (l2009tl developed 
similar models for solar-mass stars. They find a coupling 
timescale of 110 Myr for slow rotators and 6 Myr for fast 
rotators. In comparison, our models yield 30 and 10 Myr, re¬ 
spectively. While the values for slow rotators broadly agree, 
the coupling timescale for fast rotators differs widely be¬ 
tween our models. A major difference between former and 
current models lies i n the adopted wind brak ing la w. While 
iBouvie rj (l2008h and llrwin fc Bouvied (l2009t l used iKawalerl 
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Fig. 7. The angular momentum transport rate between the 
radiative core and the convective envelope is shown as a 
function of time for the three rotator models: slow (red), 
median (green), and fast (blue), and for the three mass bins 
explored here: 1 M 0 (solid), 0.8 M 0 (dotted), and 0.5 M 0 
(dashed). 


ljl988l ) braking law, we use here [Matt et ~aI1 (|2012aT ) pre¬ 
scription. As the latter braking law is more effective at ex- 
tracting angular momen tum than the former (cf. Fig.8 of 
iGallet fe Bouvieil 120 1 3D . a shorter core-envelope coupling 
timescale compensates for the enhanced surface braking 
torque and produces similar spin-down timescales on the 
MS. This instance clearly illustrates how the absolute val¬ 
ues of T c _ e depend on the assumed wind braking torque. 

The coupling timescale r c _ e we derived from our model 
also significantly increases towards lo wer mas s es, a result 
that is consistent with t hose of llrwin et al.1 ([2007) and 
iDenissenkov et all (|2010D . Figure [7] displays the evolution 
of the transport rate A J/r c _ e as a function of time for the 
three rotator models over the three mass bins considered 
here. The angular momentum transport rate of a 0.5 M 0 
star is 2 orders of magnitudes smaller than for a 1 M 0 star, 
as the coupling timescale shortens from ss 300 Myr to ~ 30 
Myr for 0.5 and 1 M 0 , respectively. 

Another way to estimate the angular momentum ex¬ 
change between the core and the envelope is to estimate 
the viscosity v i nvolved in the angu l ar mo mentum trans¬ 
port mechanism. IDenissenkov et all (|2010T > converted this 
viscosity to a characteristic timescale that can be compared 
to the core-envelope coupling timescales used here. This 
conversion can be expressed as r c _ e = 6.3 x 10 20 ^ _1 - 208 , 
with v in cm 2 s _1 . The characteristic timescale associated 
with angular momentum transport in t he st ellar i nteri or is 
related to the turbulent viscosity v dHeeer et al.ll2000i ) by 

T c - e = l 2 lv (13) 

where l is the characteristic length of the redistribution 
current. Hence, identifying the constant 6.3 x 10 20 above to 
/ 2 , we find a scale of l = 2.51 x 10 10 cm i.e. 40% i?* for a 1 
Mq star. This suggests that the diffusive approach and the 
two-zone models are somewhat equivalent and yield similar 
results. 

New constraints on the angular momentum transport 
rate in stellar interiors have emerged from the recent as- 
terosismic results derived from the analysis of Kepler light 
curves. In particular, the detection of mixed gravity and 
pressure modes have allowed various groups to derive the 
internal rotation profile of red giant stars, evolving off 



Fig. 8. Coupling timescale as a function of surface angular 
velocity at ZAMS for three mass bins: 1 M 0 (solid), 0.8 M 0 
(dotted), and 0.5 M 0 (dashed). The values derived from our 
models are shown as black dots (cf. Tabic [2]) and the red 
lines represent a power-law fit to the data. 


the main sequence fe.g. iDeheuvels et al.l 120121: Mosser 


et al. 2012f ). Using a diffusive approach. lEggenberger et al.l 
(|2012bh showed that an enhanced ’’anomalous” viscosity 
of v = 3 x 10 4 cm 2 s _1 is required to account for the 
rotational splitting of the red giant KIC 8366239, whose 
rotational profile is much shallower than diffusive models 
would predict. Using the conversion we derived above from 
IDenissenkov et abl (2010), this translates into a coupling 
timescale of 78 Myr, i.e. not very different than the one 
we derived from our PMS-MS models. This opens the in¬ 
triguing possibility that we might be dealing with the same 
physical process being instrumental from the PMS to the 
MS and beyond redistributing angular momentum in stel¬ 
lar interiors, thus reducing the amount of rotational shear 
to be expected at any evolutionary phase. 

Ideally, one would hope to derive coupling timescales 
directly from the theory of angular momentum transport 
processes. In practice, a characteristic timescale for angular 
momentum redistribution is difficult to derive from analyti¬ 
cal calculations and numerical simulations. In hy drod y nam- 
ical s imulations, such as the STAREVOL model (ISiess et al.l 
|2000D . the angular momentum redistribution timescale can 
be linked to physical mechanisms such as meridional cir¬ 
culation and thermal diffusivity. For these processes, the 
coupling timescale varies from few 100 Myr to few Gyr (L. 
Amard, private communication). For internal gravity waves 
the associated coupling timescale is about 10 Mvr (Zahar 
et al. mW- - - 

A recent attempt to deri ve coupling t i mesc ale analyti¬ 
cally has been performed bv lQglethorpe fe Garaudl (|2013f h 
Assuming that a large fraction of the radiative core is 
held in solid-body rotation by the internal magnetic field 
that is confined below the tachopause (i.e. the base of 
the tachocline), while angular momentum is transferred 
through the tachocline by large-scale meridional flows, 
they derive a coupling timescale proportional to the local 
Eddington-Sweet timescale across the tachocline finding 


T c -e — tES x 


Ip: (I core "t" Ic. on v ) 


(14) 
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with 

t ES = N 2 J2n 2 conv x (S/R core ) 2 x S 2 / Ktc , (15) 

where J tc , Kto N tc , and S are the moment of inertia, ther¬ 
mal diffusivity, buoyancy frequency, and thickness of the 
tachocline, respectively. 

With these analytic expressions, and using the param¬ 
eter values of their reference model, the average coupling 
timescale would be f c _ e = 0.34 Myr for the 1 M 0 stars. 
For the 0.8 and 0.5 M 0 , we used the values of the density, 
buoyancy frequency, and thermal diff usivity extracted from 
the STAREVOL model (|Siess et al.ll2000l . Louis Amard, 
private communication}_and^xtra£olate^ i Oglethorpe & 
Garaud (l2013f ) model to derive r c _ e =8 Myr and r c _ e =60 
Myr, for 0.8 and 0.5 M 0 , respectively, adopting the surface 
velocities of median rotators at ZAMS. We caution, how¬ 
ever, that we computed these values assuming the same 
tachocline thickness for lower mass stars as for solar-mass 
stars, while thicker tachoclines would lead to much longer 
coupling timescales (cf. Eq. 15). Even though these values 
do not match those we empirically derive from paramet¬ 
ric models, we note that theoretically predicted coupling 
timescales increase towards lower masses as required by the 
ang ular momentum evolut ion models. 

Oglethorpe & Garaud (2013j)’s derivation further pre¬ 
dicts coupling timescales that are inversely proportional 
to the square of the surface velocity, i.e. r c _ e oc l/Q 2 onv . 
Figure[8]shows the variation of the coupling timescales upon 
angular velocity as derived from our parametric models. A 
power-la w fit yields T r -e oc °' 6 ^ while a previous 

study bv lMoraux et al.l (l2013f ) of PMS and ZAMS clusters 
suggested r c _ e oc 1/fi. While these parametric slopes are 
at variance with the analytical predictions, they neverthe¬ 
less agree with the expected trend for a shorter coupling 
timescale in faster rotators. 


5.3. The mass-dependent efficiency of magnetic braking 


As sh own in Sect. 01 and discussed in ICallet fc Bouvierl 
(|2013l h the adopted braking law is a crucial parameter of 
angular momentum evolution models, as it dictates most 
of the rotational evolution on the main sequence. In the 
model presented here, the recent re sults of the numerical 
simulations from lMatt et al.l (l2012af ) and lGranmer fc Saail 
(20111) have been combined to infer the main properties of 
the surface wind generated by low-mass stars. Th i s new 
braking law contrasts with the well-known iKawalef (ll988l f 
prescription used in alm ost all angular momentum evolu- 
tion models so far (e.g. Bouv ier j [20081; Ewin &; Bouvierl 


2009: iDenissenkov et al.l 1201 ft Soada et al.l 2011 1 as well 

as with the KswalernnodifiedjnodeLjaroEOse^Jo^Reiners 
& Mohantv d2012f ). A detailed comparison between these 
brakin g laws is provided in Sect . 5.2 of lGallet, fc Bouvierl 
(|2013l f (se e also I Matt, et aT1l2014fl . 

In the iMatt et "all (2012a) prescription adopted here, 
together wi th the dynamo re l ation ship and mass-loss esti¬ 
mates from lCranmer fe Saaii (|201 11 1 . the predicted braking 
torque already includes the mass dependency. This is seen 
in Eqs. ms where mass-dependent parameters appear ex¬ 
plicitly (namely: R*, v esc , T e //, Ro , L*, p*). Hence the 
braking law takes the form 


T wind = I<1 x /(fi*, M*, JS*, .B*, M*,...), 


(16) 


where K\, the calibration factor, is expected t o be a con¬ 
stant for a given magnetic topology. Thus, in IMatt et al.l 
(l2012af ). who assumed a large-scale dipolar field, ATi=l-30. 

We were unable to fit the observed rotational evolution 
of low-mass stars with a constant wind braking calibration 
factor over the mass range explored here. As summarised 
in Table 01 the K\ parameter had to be recalibrated to a 
value of 1.7, 3.0, and 8.5, for 1 M 0 , 0.8 M 0 , and 0.5 M 0 
models, respectively. The K\ factor was, however, kept the 
same at a given mass for the slow, median, and fast rota¬ 
tor models. The increasing K\ factor towards decreasing 
mass suggest that the efficiency of the magnetic braking by 
stellar winds increases towards lower massed- We offer two 
possible explanations for this result. 

The stellar wind simulations of lMatt et al.l ( 2012a l com¬ 
puted for a dipolar magnetic geometry were expanded to 
qua drupolar and octupolar magnetic topologies bv Reville 
et al. (|2015f i. As expected, all other parameters being equal, 
the radially averaged Alfven radius is found to be sys¬ 
tematically smaller for more complex magnetic geometries. 
Hence, our empirical recalibration of the K\ factor could 
be interpreted as indicating that higher mass stars have 
more complex magnetic geometries than lower mass stars, 
thus reducing the efficiency of magnetic braking by stel¬ 
lar winds at higher masses. A mass-dependent magnetic 
topology is indeed suggested for main-sequence stars, with 
solar-mass stars having more multipolar magnetic structure 
than lower mass star s whose magnet i c top olog y tends t o be 
dipole dominated dPetit et al.1 120081: iDonati fe Landstreetl 


120091 : iMorin et ahll201fll : iGretrorv et alJ 2012f l. 

A more prosaic explanation, however, may simply be 
that the mass-loss rate prescription that we use depends 
too strongly on stellar mass. As seen from Eq. 0] and 0] we 
find 


r wind ex K\ x MItT ^ K l x (17) 

so that we could produce the same results by varying 
the mass-loss rate instead of recalibr atin g the K\ facto r. 
Keeping the K\ factor to 1.30 as in IMatt, et ~aT1 (l2012alf 
would require multiplying the mass-loss rates of 1 M 0 , 0.8 
M 0 , and 0.5 M 0 stars by factors of 2.6, 20, and 860, respec¬ 
tively. Compared to the mass-loss rates used in the para¬ 
metric models above (see Fig.QJ, which range from 10 -16 5 
M 0 yr _1 for 0.5 M 0 stars to 10 -13 Mf^yr” 1 for young solar- 
mass stars, this would translate into mass-loss rates in the 
range ~10 _14 -10 -13 M 0 yr _1 over the whole mass range. 
The latter values might actuall y be in b etter a greement 
with the observations (e.g. IVidotto et al.ll2013 l). Clearly, 
additional measurements of magnetic topologies and mass- 
loss rates are needed for young solar-type stars to be able 
to better calibrate the braking torque they undergo as they 
evolve on the early main sequence. 

As discussed in Sect. ro lower mass stars take longer 
to spin-down on the main sequence, but once they have 
reached rotational convergence, their spin rate is lower than 
that of solar-mass stars (cf. Fig. [5]). The longer spin-down 
timescale and lower terminal velocities of lower mass stars 
both directly results from the braking law. Since stars in 
the 0.5-1.0 M 0 mass range have about the same maximum 
velocities at ZAMS, the spin-down timescale primarily de¬ 
pends on how much angular momentum is extracted in the 

5 Note that the braking torque always remains higher for 

higher mass stars, as shown in Fig. [4] 
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Fig. 9. The spin-down rate Cl = J/1* is plotted as a func¬ 
tion of angular velocity for 1 M 0 (solid), 0.8 M 0 (dotted), 
and 0.5 M 0 (dashed) stars. Note the cross-over of the de¬ 
celeration curves as the stars go from the magnetically sat¬ 
urated to the unsaturated regimes. 


magnetically saturated regime. As shown in Fig. 3J this 
quantity decreases towards lower mass stars, thus account¬ 
ing for their longer spin-down timescale. The terminal ve¬ 
locity, however, depends mostly on the extraction of angular 
momentum once the star is nearing rotational convergence, 
i.e. in the unsaturated regime (low rotation rate). Fig. [I] 
shows that the braking torque is then about the same for 
all stars in the 0.5-1.0 M 0 mass range, and since the lower 
mass stars have a lower moment of inertia, they are more 
easily braked, eventually reaching lower terminal velocities. 
Hence, the two major properties of the rotational evolution 
of low-mass stars on the main sequence are fully encoded 
into the mass-dependent braking law. This is illustrated in 
Fig. [9] which shows the deceleration rate of low-mass stars, 
= J/I, as a function of their instantaneous velocity. It is 
seen that the curves for the three mass bins cross-over at in¬ 
termediate velocity, i.e. higher mass stars are more rapidly 
spun down at high velocities, while lower mass stars are at 
low velocities. 

5.4. Gyrochronology and the Skumanich relationship 

Figure [T] shows the observed period-mass distributions of 
the various clusters used in this study together with the 
model predictions for evolution of slow, median, and fast 
rotators over the 0.5-1.0 M 0 mass range. It is seen that 
the stellar spin smoothly evolves with time at all masses 
investigated here and models closely follow this evolution. 
Hence, the possibility of uniquely relating the spin rate to 
the age of the star i.e. to explore gyrochronology as orig¬ 
inally sug gested by iBarnesI (2003 1 following the pioneering 
work of I S kumanich! (Il972f ). The Skumanich’s relationship, 
fl* oc t -1 / 2 , was derived for solar-mass stars with an age 
between 100 Myr and the Sun’s. In the models presented 
here, the asymptotic behaviour of low-mass stars at large 
ages appear to follow this relationship quite closely (cf. Fig. 
[5]). It is however important to stress that: i) this relation¬ 
ship is valid only past the time of rotational convergence, 
i.e. after ^0.5 Gyr for solar-mass stars and after ~3 Gyr 
for 0.5 M 0 stars; ii) prior to the epoch of rotational conver¬ 
gence, there is no one-to-one relationship between a star’s 


age and its surface velocity as the initial dispersion of ro¬ 
tation rates has not been erased yet; and iii) the rotational 
evolution on the early MS can be steeper (fast rotators) or 
shallower (slow rotators) than the Skumanich relationship 
would predict depending on the amount of angular momen¬ 
tum being stored in the stellar core and slowly resurfacing 
into the envelope (cf. Fig. [5]). 

On the late-MS, all stars in the 0.5-1.0 M 0 mass range 
have reache d rotati onal convergence. Indeed, the recent re¬ 
sults of iMeibom et al.l (12015! ) on NGC 6819, a 2.5 Gyr 
open cluster, yield a tight mass-rotation period relation¬ 
ship over the mass range 0.85-1.3 M 0 at this age. In this 
cluster, solar-mass stars have rotational periods narrowly 
distributed within PiM 0 =17-19d; while 0.85 M 0 stars 
have periods in the range Po.85M 0 =21-24d. Our solar-mass 
model interpolated at this age predicts F > ijvr 0 =18.5±0.8d 
(with a median period of 17.7d), while the 0.8 M 0 model 
yields Po.8Af 0 =24.8±2.Od (median P=22.67d), in excellent 
agreement with observations. 

The larg est set of rotational per iods has been recently 
obtained by iMcQuillan et al.l (2014! ) for a sample of more 
than 34000 field stars using the Kepler satellite. We repro¬ 
duce their results in Fig, fill where the spin rates are plotted 
as a function of mass, and where we over plotted isochrones 
predicted from our 0.5, 0.8, and 1.0 M 0 models. It is seen 
that the lower envelope of the Kepler rotational distribu¬ 
tion is well reproduced by the models at an age of 7 Gyr, 
as expected for the oldest field stars of the Kepler sample. 
The upper envelope of the main rotational locus roughly 
corresponds to 0.5 Gyr-old stars, and the ’’rain” of stars 
above it consists of even younger field dwarfs. 

Finally, we should emphasise that there is no dichotomy 
in the rotational evolution of low-mass stars. We illustrated 
here models for slow, median, and fast rotators and all three 
models rely on the same assumptions. Indeed, the evolution 
of slow and fast rotators alike is driven by the same under¬ 
lying physical processes. These processes act on different 
timescales depending on stellar mass and initial velocity, 
which leads to the observed, peculiar and rapidly evolving, 
shape of the cluster’s rotational distributions (cf. Fig. !T()1) 
and, ultimately, to the field distribution (cf. Fig. [IT]). Thus, 
the observed rotational distributions of young open clus¬ 
ters and field stars do not prompt for different processes 
operating in fast and slow rotators. Instead, the sequence 
of rotational distributions from 1 Myr to several Gyr can be 
described as deriving from the smooth temporal evolution 
of a single, widely dispersed initial distribution of angular 
momenta at birth. 

There are other s two-z o ne models in the lit e rature , such 
as Keggense^alJ (1995! ) . Ilrwin fc Bouvier] (20091) . and 


iDenissenkovI ( 2ofol ) . These models are based on the lKawalei! 
braking law and their numerical structures are quite 


similar to our model (star-disk interaction, core-envelope 
decoupling, and wind braking). The mai n differe nce lies 
in the wind braking law used. The iMatt et al.1 (2012a!) 
braking law contains more physics than Kawaler’s, espe- 
cially thr ough the mass-loss rate, which is not included in 
[Kawalerl (|1988f ). Moreover, while in these models the sat¬ 
uration of the braking low has to be manually set, in the 
model presented in this article the saturation is fixed by 
the observations via the magnetic filling factor. 

We also highlight that there are other modelling ap¬ 
proaches that do not necessarily invoke the two-zone for¬ 
malism. These models usually assume solid body rotation 
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Fig. 10. Same as Fig. [T| The red and blue crosses are the values derived from our evolution model for the slow and fast 
rotator models, respectively. 


for the whole stars. [Barnes! (I2010T ) and lBarnes fe Kiml (I2010T) 
propo se a model based on gyrochronology analysis (I Barnes! 
[200311 and a Kawaler-modified braking law. This model 
only works for moderate-to-old MS stars once the conver¬ 
gence is achieved, but fails to reprodu ce the ro tation be¬ 
haviour of young stars. In [Reiners fe Mohantvl (120121 ) the 
braking torque is inversely proportional to the mass-loss 
rate, which is assumed to be constant during all the evo¬ 
lution. This behaviour somehow contrasts with what we 
expect, i.e. that the magnetic braking is induced by the 
quantity of angular momentum extracted by the stellar 
wind. However, their model can easily be extended to very 
low-mass sta r th anks to th eir mass and radius dependen¬ 
cies. Fi nally. Brownl (2014) proposed a more sophisticated 
[Barnes! ((20101 ) like model. His models are quite good for 
intermediate ages (~ 200 Myr) but fail to reproduce the 
youngest clusters (a Per) because of a lack of slow rotation 
(due to the solid body rotation) and the old cluster (M37, 
NGC6811) especially for low-mass stars. 


6. Conclusion 

The rotational evolution of low-mass stars in the range from 
0.5 to 1.1 Mq can be described, from birth to the end of 
the main sequence by parametric models that rely on a 
limited number of physical processes: star-disk interaction, 
wind braking, and core-envelope decoupling. The physical 
processes involved in the angular momentum evolution are 
included in the models through parametric prescriptions 
that either rely on observational evidence (e.g. rotational 
regulation during PMS star-disk interaction), are physically 
driven simplifications of actual processes (e.g. core-envelope 
decoupling), or are based on recent numerical simulations 
(e.g. wind braking). We explored here how the model pa¬ 
rameters vary with mass to get deeper insight into the un¬ 
derlying physical mechanisms. The disk lifetimes and initial 
rotation periods seem to have little dependence on mass, 
and the former are consistent with the distribution of disk 
lifetimes derived from the evolution of IR excess in young 
stars. However, we do find a correlation between initial pe¬ 
riod and disk lifetime, which may point to the impact of 
protostellar disks in establishing the initial distribution of 
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Fig. 11. Angular velocity (left axis) and rotation periods (rig ht axis ) as a function of mass for the 34030 field stars whose 
rotational pe riods have be en mea sured bv iMcQuillan et alJ (120141 ). The stellar mass was derived from T e // using the 
isochrones of lBaraffe et al.l (119981 ) for a mean age of 1 Gyr. 


angular momenta. In contrast, the core-envelope coupling 
timescale and the wind braking efficiency strongly vary with 
mass, and both increase towards lower mass stars. While 
the former result may reflect the properties of a still to 
be identified angular momentum transport process oper¬ 
ating in stellar interiors, which is also shown to strongly 
depend on rotation, the latter may derive from a change 
in the magnetic topology of dynamo fields as one goes to 
lower masses. We find that these models do reproduce the 
run of rotational period distributions as a function of time, 
from the youngest star-forming regions to the oldest open 
clusters, and provide a remarkable fit to the oldest field 
stars from the Kepler sample. The implications are mani¬ 
fold. A large amount of angular momentum must be stored 
in the radiative core during several 100 Myr on the early 
main sequence, which will undoubtedly impact on the stel¬ 
lar properties, such as lithium content, even long after the 
end of the spin-down phase on the main sequence. Also, 
the build-up of a wide dispersion of rotational velocities 
at ZAMS from an initially dispersed PMS distribution and 
its subsequent evolution on the early MS partly reflect this 
process. The models also naturally account for the longer 
spin-down timescale of lower mass stars and their lower 
velocities compared to solar-type stars at the end of the 
spin-down phase. Finally, we show that the models nat¬ 


urally yield a Skumanich-type rotational evolution on the 
late main sequence for stars in this mass range, but that un¬ 
til about 1 Gyr, there is no one-to-one relationship between 
stellar rotation and stellar age, thus undermining the deter¬ 
mination of accurate stellar ages for individual young stars 
from their rotational properties. Indeed, the slow release on 
the early main sequence of angular momentum hidden in 
the stellar core delays the epoch at which a Skumanich-type 
rotational evolution eventually sets in. 
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